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Abstract 

The 1 /N expansion solutions for the interacting boson model are extended 
to higher orders using computer algebra. The analytic results are compared 
with those obtained from an exact diagonalization of the Hamiltonian and are 
shown to be very accurate. The extended formulas for level energies and E2 
transitions will be useful in the analysis of high-spin states in both normal 
and super deformed nuclei. 
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Application of the interacting boson model |Q (IBM) to the high-spin states of 
normal and superdeformed nuclei has been rather slow, and certainly not commensu¬ 
rate with the current experimental activity in this held. The reasons for this are i) the 
necessity of including g bosons in the basis in addition to the usual s and d bosons, 
and ii) the need to include broken pairs (quasi-particles) in the basis if backbending 
occurs. In this work, we address the former problem and refer to Ref. |0] for recent 
developments in the latter. The importance of g bosons has been well established 
in studies involving low-lying levels of deformed nuclei (see Refs. ii for reviews). 
There are also many experiments on high-spin states where results are compared 
to the sd-IBM calculations showing its inadequacy. The difficulty with performing 
sd^f-IBM calculations for deformed nuclei (with boson numbers iV = 12 — 16) is that 
the basis space is too large, and numerical diagonalization of a Hamiltonian is not 
possible even on a supercomputer. Therefore various truncation schemes have been 
devised. These schemes are geared towards reproducing the low-lying spectrum and 
can not be expected to give reliable results for the high-spin states. The situation 
gets even worse in the case of superdeformed nuclei where recent microscopic studies 
have shown that the number of bosons is around iV = 30 — 40 p, Q , which renders 
any attempt to numerical diagonalization futile. 

A more promising approach is to use the angular momentu”m projected mean 
held theory which leads to a 1/A^ expansion for all matrix elements (m.e.) of interest 
[0. Initial calculations in the 1/N expansion formalism were restricted to order 
1/N^ which was sufficient for low-spin {L < N) phenomenology. For an accurate 
description of high-spin properties (e.g. dynamic moments of inertia), however, 
one needs to include terms up to order 1/iV® (note that these terms actually go 
as (L/iV)®, and are not small for high-spin states). Hand calculation of m.e. to 
such high orders is impractical and has never been attempted. Recent advances 
in computer algebra have hnally broken this impasse. The purpose of this letter 
is to report on the results of such a 1/A^ expansion calculation facilitated by the 
Mathematica software [§. Applications are made to the superdeformed bands in 
Hg isotopes. 

We consider a general formulation of the IBM and introduce the boson creation 
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and annihilation operators bl^,bi^ with I = 0,2,4,.... The gronnd band can be 
written as a condensate of intrinsic bosons as 


|]V,x) = (iV!)-‘/2(),t)A'|o), 6* = 


( 1 ) 


where xi are the normalized boson mean helds which are associated with the de¬ 
formation parameters of the system. Here we assnme that the system is axially 
symmetric, which is a good approximation as will be seen later. The other bands 
can be obtained from the gronnd band by acting with the orthogonal intrinsic boson 
operators. A general boson Hamiltonian with one- and two-body terms is given by 


2ln 


H = Y.eihi+ Y. 

I k=0 

M jl 


( 2 ) 


where brackets denote tensor conpling of the boson operators, bi^ = (— 
and hi and are the boson nnmber and mnltipole operators, respectively. The 

parameters in the model are the single boson energies the mnltipole strengths 
Kfc, and the coefficients tkji- For consistency, the same mnltipole operators are nsed 
in the calcnlation of electromagnetic transition rates. 

Matrix elements of any operator with angnlar momentnm projection can be 
evalnated efficiently nsing boson calcnlns and angnlar momentnm algebra techniqnes 
[0. We refer to the original references for details and give here the intermediate 
resnlts that form the basis of the Mathematica calcnlation. The exact expectation 
valnes of the nnmber and mnltipole operators (@) in the gronnd band are given by 

{Ai)l = J2(L0100\I0)^F(N -1,1), 

Ijm. TtF), = AiA! j2{Lom\io)^F(N -1, /) 


X -I 21 + 1 J 

+N{N-1) Y tkjitkjn'XjXiXfXr{j0f0\J0){M'0\J0) 
jij'i'j 

J ' 


X 


J J 
V I 
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Y{L0J0\I0)^F{N -2,1) 


(3) 


Here F{N,L) denotes the normalization integral which has the 1/N expansion P] 

.T'", (4) 
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where bar denotes the angular momentum eigenvalues, L = L(L +1) and a = 
represents the “average angular momentum squared” carried by a single boson. The 
coefficients anm in Eq. have been calculated recently utilizing the Mathematica 
software |0 . It is clear from the expressions that, in principle, one can evaluate 
these m.e. to any desired order in 1/A^ (see Ref. |^). However, the amount of 
algebraic manipulation required grows exponentially with each order, and one is 
soon forced to give up. In contrast, Mathematica is very efficient at such “algebra 
crunching”. 

Before presenting the hnal results, it will be useful to comment on the general 
form of the m.e. of a /c-body operator O, and illustrate the concept of layers in the 
1/N expansion (see Ref. 0 for more details) 


{o)l=n'^y: 


Or, 


The expansion coefficients Onm in Eq. (|^) involve various quadratic forms of the 
mean helds xi corresponding to the single-boson m.e. of O and its moments. The 
terms with n + m — 1 = i, constant, are referred as the ith layer. An expansion in 
layers rather than in 1/N is preferred on both practical and physical grounds. This 
is because terms in the same layer have similar forms (although they have different N 
dependence), while complexity of terms grows exponentially with increasing layers. 
For example, a complete calculation to order 1/N^ would involve terms belonging 
to the fourth, hfth and sixth layers which are very complicated yet completely un¬ 
necessary as their contribution to the m.e. is beyond the desired level of accuracy. 
Here we restrict ourselves to a third layer calculation of the m.e. which is found to 
be sufficiently accurate for description of high-spin states. For the one-body m.e. in 
Eq. (!]), one obtains 


-V 


(5) 


(ffi)L — Nxf^l + + (1 “ ai/a)l + F/2j 

-|- - -^ (cb 2,(1 — 7cii/3 — (Kii -|- 5oi /4o — 02 /3 

[aNy^ ^ 

-|-(—1 — 2(1 -|- 2(ii ^(X\j2(x — 5ci^/2ci^ -|- (X‘i/2(xyi 

-|-(—7/6 — a + 5ai/4a)P — F/6^ 


L 
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1 


(2a + 2 oF 


2ai + (—2 — 2a -I- ?>ai/a)l — 
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3a — 12a^ — 4a^ + 21ai/2 + llaoi — 15a?/2a + 802/2 

[aNy ^ 

+ (3 + 12a + 4a^ — 33ai/2 — 14ai/a + 25a^/2a^ — 2a2/a)l 
+ (7/2 + lla/2 - 5ai/a)P + P/2)' 


2(aiV)^ 


aN 


L3 


3(aiV)f 


(—a — 2a + 3ai/2 + (1 + 2a — 2ai/a)l + 1/2) 

— (^4a + 21a^ + 14a^ — 16ai — 51aai/2 + 13a^/a — 2a2 

+ (—4 — 21a — 14a^ + 34ai + 20ai/a — 390^/20^ + 5a2/2a)l 
+ (-4 - 17a/2 + 13ai/2a)P - P/2)' 

—a — 6a^ — 6a^ + 25ai/6 + 9aai — 15a^/4a + 5a2/12 
+ (1 + 6a + 6a^ — 45ai/4 — 5ai/a + 21a^/4a^ — 02/20)! 

+ (5/6+ 9a/4-3ai/2a)P +P/12]}, (6) 


where a„ = denotes the higher moments of a. We have checked Eq. 

(H) against two resnlts: i) it satishes the nnmber conservation, i.e. = N, 

ii) it reprodnces the analytic formnlas available in the SU(3) limit [|l|. A similar 
Mathematica evalnation of the two-body m.e. in Eq. (S) gives 


(TW . = N^[Uk + - Ukl + aCk) 


(aA^)2 
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{aNf 


}(—2a + ai)Uk + (1 — a — ai/a)Uki + Uk 2 / 2 , + a/Ck — ciCki^ 
(^{2a + 2a^ — 14ai/3 — aoi + 5a^/2a — 2a2/3)Uk 


+ (—1 + a — a\/2 + 7ai/2a — 5a]^/2a + a2/2a)U) 


kl 


{aNY^ 


+ (-7/6 + 5ai/4a)f/fc2 - Uks/Q 
+ (—a^ + aai/2)Ck + (a — ai)^^! + aCk 2 / 2 ^ 

-2aUk + Uki 


1 


aN 
1 


H—rr ((4a + 2a^ — 4ai)f/fc + (—2 + a + 3ai/ a)Uki — Uk 2 — oYC^ + oCfci) 
:((—6a — 16a^ — 4a^ + 21ai + 15aai — lba\/a + 3a2)Uk 
+ (3 + 2a — 2a^ — 4ai — 14ai/a + 2ba\/2aY — 2a2/a)Uki 
+ (7/2 + 2a — 5ai/ a)Uk 2 + Uk‘i/2 


{aNY 
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+ (2a^ + 2a^ — 2aai)Ck + (—2a — 2a^ + 3ai)Cfci — aC, 


k2 


2{aNy 




3(aiV)f 


(—2a — 2a^ + 3ai)Uk + (1 — 2ai/a)Uki + f/fc 2/2 

-f(8a + 30a^ + 14a^ — 32ai — 37aai + 26a?/a — 4a2)f/fc 

aN ^ 

+ (—4 — 8a + 2a^ + 29ai/2 + 20ai/a — 39ai/2a^ + 5a2/2a)Uki 
+ (—4 — 9a/2 + 13ai/2a)t/fc2 — Uk3/2 
+ (—a^ — 2a^ + 3aai/2)C'fc + (a + 2a^ — 2ai)Cki + aC*fc 2 / 2 ^ 

(—8a — 36a^ — 24a^ + lOOai/3 + 54aai — 30a^/a + 10a2/3)Uk 
+ (4 + 12a — 24ai — 20ai/a + 21ai/a^ — 2a2/a)Uki 


+ (10/3 + 6a — 6ai/ a)Uk 2 + Uks/S 


(7) 


Here the quadratic forms Ckn arise from normal ordering and simulate an effective 
one-body term 

Ckn = {2k + 1) ^ C{tkjixif/{21 + 1), (8) 

ji 

while Ukn represent the genuine two-boson interaction 
Ukn — 

S'l 7''(j0j'0|/0)(/0/'0|/0) I ^j^^tkjitkj'i'XjXiXfXi>. (9) 

For a given multipole, these sums can be evaluated in closed form using Mathemat- 
ica. For the quadrupole interaction, the hrst four terms needed in Eq. (0) are given 
by 


U 2 = Ay U 21 = (2Hi - 3A)A, 

U 22 = (2H.2 — 24H.1 + 18H.)H + (Hu — H 2 + 7Hi)Hi + (Hu — H2)^/12, 
U 23 = (2H3 - 36H2 - I8H11 + 192Hi - 144H)H 

+ (3H2i — 3 H 3 + 48 H 2 + 24Hii — 146Hi)Hi/2 
— ( 3 H 21 — 3 H 3 + 25 H 2 — 14Hii)H2/12 


+ (3H2i — 3H3 + llHii)Hii/12, 


where the quadratic forms A^n in (P^ are dehned as 


Amn = ^ {j0l0\20)t2jlXjXl, 
jl 


( 10 ) 


( 11 ) 
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and correspond to various moments of the single-boson quadrupole m.e. (note that 
the zero subscripts are suppressed for convenience). The quadrupole m.e. given by 
TO reproduces the well known Casimir eigenvalues in the SU(3) limit, hence also 
passes the SU(3) test. 

To obtain the hnal energies, one has to perform variation after projection (YAP) 
and determine the mean helds xi for each L. This can be done algebraically using an 
ansatz for xi similar to Eq. (|^), which leads to sets of linear equations for the higher 
order terms |0|. As the algebraic solution of the VAP equations are rather lengthy, 
we defer their discussion to a longer paper. Alternatively, one can determine xi 
numerically using, for example, the simplex method. Both methods lead to equally 
accurate results. Expressions for the excited band energies are obtained in a similar 
manner, though they are somewhat more complicated due to contributions from 
orthogonality and band mixing effects, and will be given in a future paper. 

The other experimentally relevant quantities in the study of high-spin states are 
the E2 transitions. In this case, the hrst layer results are found to be sufficiently 
accurate, hence higher order calculations are not necessary. For completeness, we 


quote the E2 transition m.e. in the ground band |10 


{E II T(2) II L) = NL{L0 20|L'0){A - —(Ai - 3A - aA) 

1 


S{aNy 


(Aa - All - lOAi -F 12A -F 4aA)(L' -F L) 


( 12 ) 


~g(A 2 — All ~ 6 Ai -|- 6 A)(L' — L)^ |, 

where L = [2L -|- 1]^/^ and Amn are dehned in (pA)). Further expressions for in- and 
inter-band transitions among the ground, (3 and 7 bands, including band mixing 
effects, are given in Ref. [O. A fortran code for the 1/N expansion calculation of 


energies and E2 transitions is available from the authors upon request jll 


Before applying the 1/N expansion results, we compare them with those obtained 
from an exact diagonalization of the Hamiltonian |^. For this purpose we choose 
two different parametrizations of sd^f-IBM: (a) a pure quadrupole interaction which 
is appropriate for superdeformed states [ 0 , |^, (b) quadrupole interaction plus a g 
boson energy which reflects weaker coupling of g bosons in normal deformed states. 
In both cases, the quadrupole parameters g 22 ,? 24 ,Q '44 are scaled from their SU(3) 
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values with a single factor q (go 2 = 1, Qji = hji)- The q factor gives a simple 
measure for the breaking of the SU(3) symmetry which is necessary for a realistic 
description of deformed nuclei. Of necessity, the boson number is hxed at iV = 10, 
which requires some scaling of the parameters from their realistic values at larger 
N. Fig. 1 presents comparisons for the ground band energies. The top hgures 
show El/L as a function of L = L{L + 1). In case (a), the agreement between the 
third layer 1/N expansion calculation (solid line) and the exact results is excellent 
up to the maximum spin L = AN. The second layer results (dotted line) start 
deviating from the exact ones around L = 2N as do the third layer results with 
VBP (dashed line). In case (b), the agreement is still very good, with only a few 
percent difference at the very high-spin region. Here, the VBP calculation shows 
even larger deviations, underscoring the importance of the YAP procedure for high- 
spin states. In the bottom hgures, we compare the dynamic moment of inertia 
which is much more sensitive to changes in structure. Thus, the inadequacy of 
the second layer calculations, which is not very clear in the energy plots, becomes 
obvious in plots. In Fig. 2, we compare the hrst layer results for the E2 
transition m.e. ([I2|) with the exact ones. The agreement is again very good up 
to very high-spins in both cases. The test cases discussed above indicate that the 
extended formalism can be applied with conhdence in the spin region L = N-3N 
which covers the presently available high-spin data. This, incidentally, also shows 
the adequacy of the axial symmetry assumption made in the beginning. 

The analytic expressions derived for energies and E2 transitions will be useful 
in the study of high-spin states in both normal and superdeformed nuclei. Here, we 
present an application of the 1/N expansion formalism to superdeformation which 
is more topical and harder to treat by numerical diagonalization. In super IBM, 
as proposed by Otsuka and Honma [^, ||, normal bosons are supplemented with 
superdeformed bosons which correspond to the Cooper-pairs in the superdeformed 
potential. The number of superdeformed bosons, Nguper is typically around 30-40, 
and because of large deformation, g boson effects are important. Thus, the super 
IBM offers a fertile ground for the application of the 1/N expansion. We use the 
energy formula (|^ to describe the superdeformed bands in the Hg isotopes. The 
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dynamic moments of inertia, that resnlt from the qnadrnpole Hamiltonian are 
shown in Fig. 3. The three qnadrnpole parameters ^ 22 , ^' 24 , <?44 are scaled from their 
SU(3) valnes with a single factor q. Nsuper is determined from microscopic calcnla- 
tions, 1^, P and n and q are htted to the experimental data. A good description of 
experimental (circles) is obtained. We note that the SU(3) limit corresponds 
to a rigid rotor and wonld give a flat line for This happens becanse in the 

SU(3) limit, the mean helds xi are constant (independent of L), and the structnre 
does not change with rotation. In reality, one expects a gradnal change in dne 
to loss of pairing. The above stndy shows that this can be simulated by breaking of 
the SU(3) symmetry which results in migration of the mean helds from s to d, and 
to g bosons with increasing spin. The q values obtained in the above hts indicate 
that this breaking is around 30%. It has been suggested that the identical band 
phenomenon may be due to an underlying symmetry |T^, |T^ . It would be of inter¬ 
est to pursue this suggestion by extending the present calculations to other bands 
and also to odd nuclei. 

Another area where the 1/N expansion formalism could contribute signihcantly 
is the study of high-spin states in normal deformed nuclei. In many experiments 
on high-spin states, results were compared to the sd-IBM calculations with nega¬ 
tive connotations. This is presumably due to lack of sd^f-IBM calculations which, 
hopefully, will become more accessible with the analytic formulas presented here. 

This research was supported by the Australian Research Council, and in part 
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the Fujitsu VP of the ANU Supercomputer Facility. S.K. thanks T. Otsuka and 
M. Honma for valuable discussions on super IBM, and the members of the nuclear 
theory group at the University of Tokyo for their hospitality. 
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Figure captions 


Fig. 1. Comparison of the ground band energies (top) and dynamic moments of 
inertia (bottom) obtained from the 1/N expansion with the exact numerical results 
(circles). The different lines refer to the third layer calculation with YAP (solid 
line), third layer with VBP (dashed line), and second layer with YAP (dotted line). 
The YBP results for deviate strongly from the exact ones and are not shown 
to avoid cluttering of the hgures. The parameters of H are k = —20 keY, g = 0.7 
for the quadrupole interaction in both (a) and (b), and Sg = 500 keY in (b). 

Fig. 2. Comparison of the E2 transition m.e., Eq. (12) with the exact diagonaliza- 
tion results (circles). Solid line refers to YAP and dashed line to YBP. 


Fig. 3. Comparison of the experimental dynamic moment of inertia in 


190-194 


Hg 


(circles) with the super IBM calculations (solid lines). The data are from El 
The parameters used in the hts are Nguper = 29, 30, 31, k = —35, —34, —33 keY, 
q = 0.68, 0.72, 0.72 for 190 - 192-194 respectively. 
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